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ABSTRACT 



The Mira variable R Hydrae is well known for its declining period, which Wood & Zarro 
(1981) attributed to a possible recent thermal pulse. Here we investigate the long-term period 
evolution, covering 340 years, going back to its discovery in AD 1662. The data includes 
photometric monitoring by amateur and other astronomers over the last century, and recorded 
dates of maximum for earlier times. Wavelets are used to determine both the period and semi- 
amplitude. We show that the period decreased linearly between 1770 and 1950; since 1950 the 
period has stabilized at 385 days. The semi-amplitude is shown to closely follow the period 
evolution. Analysis of the oldest data shows that before 1770 the period was about 495 days. 
We find no evidence for an increasing period during this time as found by Wood & Zarro. 
We discuss the mass-loss history of R Hya: the IRAS data shows that the mass loss dropped 
dramatically around AD 1750. The evolution of the mass loss as function of period agrees 
with the mass-loss formalism from Vassiliadis & Wood; it is much larger than predicted by 
the Blocker law. An outer detached IRAS shell suggests that R Hya has experienced mass- 
loss interruptions before. The period evolution can be explained by two models: a thermal 
pulse occuring around AD 1600, or an non-linear instability leading to an internal relaxation 
of the stellar structure. The elapsed time between the mass-loss decline giving rise to the outer 
detached shell, and the recent event, of approximately 5000 yr suggests that only one of these 
events could be due to a thermal pulse. Further monitoring of R Hya is recommended, as both 
models make strong predictions for the future period evolution. We argue that R Hya-type 
events could provide part of the explanation for the rings seen around some AGB and post- 
AGB stars. Changes in Mira properties were already known on a cycle-to-cycle basis, and 
on the thermal-pulse time scale of ~ 10 4 yr. R Hya shows that significant evolution can also 
occur on intermediate time scales of order 10 2 -10 3 yr. 

Key words: stars: individual: R Hya - stars: AGB and post-AGB - stars: oscillations - stars: 
mass-loss - stars: variables: other - history and philosophy of astronomy 



1 INTRODUCTION 

R Hya is an unusual Mira variable. Miras are long-period variables 
found near the tip of the Asymptotic Giant Branch (AGB). They 
show mono-periodic light curves with large visual amplitudes of 
more than 2.5 mag. The periods are typically 200-500 days; the 
amplitude and shape of the light curve can vary over time but the 
periods tend to be stable. Optical data covering a century or more 
confirm the remarkable stability of the Mira pulsations (e.g. Sterken 
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et al. 1999). But in sharp contrast to this rule, the period of R Hya 
has been declining steadily for over a century^]. 

Although period jitter o f a few per cent is common among Mi- 
ras (Lombard & Koen 1993), possibly related to small changes in 
the shape of the light curves, there are only a few examples of sig- 
nificant period evolution. Other types of changes appear to be more 
common, but stars which exhibit them are automatically classified 
as semiregular (SR): the Mira classification requires pulsation sta- 
bility. The SR class is a mixture of hidden Miras and non-Mira 
stars. Examples of the former include R Dor (located on the Mira 
PL relation) which shows sud den switches betwee n a period of 
330 days and one of 180 days (Bedding et al. 1998), indicative of 



Olbers (1841 ) first noted the irregularity of the period. 
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mode switching. V Boo has shown an almost complete disappear- 
ance of its Mi ra pulsation over a cen tury, albeit without any change 
in its period ( [Szatmary et al, 1996 ). But only R Aql is known to 
show a continuous period decline similar to that of R Hya. 

Early AGB stars contain a helium-burning shell. But during 
the last 10% of the AGB, when the helium becomes exhausted, the 
shell switches to hydrogen burning, punctuated by reg ular helium 
flashes: the therma^pukesjVassfliadis^&^oo^^^^B oothroy d 
& Sackmann 1988). |Wood & Zarro (1981 ) argue that a recent ther- 
mal pulse could explain the period change of R Hya, if the star is 
presently in the luminosity decline following the peak of the pulse. 
This interpretation has generally been followed in other papers dis- 
cussing period evolution, e.g. on R C cn and T UMi ^Hawkin s, Mat- 
tei & Foster 2 001; IWhitelock 1999); Mattei & Foster 200C|; G al & 
Szat mary 1995). In support of their interpretation, W ood & Zarro 
(1981) find that the earliest observations of R Hya indicate an in- 
creasing period, which they explain with the luminosity increase 
immediately after the onset of the thermal pulse. 

In this paper we analyse data of R Hya going back to its dis- 
covery in AD 1662. The light curve is subjected to wavelet analy- 
sis, which shows how the period and amplitude (the latter available 
only since 1900) have evolved over time. We find that the decline in 
period is accompanied by a decline in amplitude. We also find that 
the period is no longer decreasing, having stabilised at 385 days in 
about 1950. The period was about 495 days before 1770; and we 
do not confirm the reported early period increase. R Hya appears 
to have evolved to its present stable period over approximately 200 
yr. 

The paper is organised as follows: In Section 2 we describe 
the data and the analysis methods. Section 3 contains a detailed 
discussion of the period evolution. In Section 4 we discuss the char- 
acteristics of the star. Section 5 discusses the pulsation evolution in 
terms of proposed AGB relations. Section 6. 1 discusses the mass- 
loss history and Section 6.2 describes the two models which can ex- 
plain R Hya-type behaviour, and show that the circumstellar rings 
observed. Finally, in Section 6.3 er discuss a possible relation to the 
rings observed around AGB and post-AGB stars. The conclusions 
are summarized in Section 7. 



2 OBSERVATIONS 
2.1 The data 

Many bright long-period variable stars have been monitored by 
amateur astronomers. The observations can be found in public 
archives^ These archives are valuable resources, and complement 
high-precision photometry datasets (e.g., Hipparcos, MACHO) that 
only cover a few years. 

The magnitudes are determined by eye, using reference fields 
that contain stars with a range of known magnitudes: a magnitude 
for the target star is established by comparison with this standard 
sequence. The accuracy is typically 0.1 mag. For red stars, sys- 
tematic differences may exist between observers. When using such 
data, sufficient observations from each individual source should be 
available to test for individual accuracy and systematic offsets. 

For R Hya (HR 5080; HIP 65835), we used data from the 



2 AFOEV and VSOLJ data are available for immediate download. The 
BAAVSS data can be requested by e-mail. The AAVSO data can be 
dow loaded for post- 1969 observations and other data can be requested by 
e-mail. 



American Association of Variable Star Observers (AAVSO), the 
Variable Star Observers League of Japan (VSOLJ), the Associ- 
ation Francaise des Observateurs d'Etoiles Variables (AFOEV) 
and the British Astronomical Association, Variable Star Section 
(BAAVSS). Only data from individual observers contributing 30 or 
more observations were used and we did not attempt to correct for 
offsets between observers. Fig. [I] shows the combined data, binned 
to 10-day averages. 

Before about AD 1890, the compilations of C annon & Pick- 
ering ( 1 9097^n ^MuTter*^Tiartwig^l^l^P ^rve"derived dates of 
maxima (and, more rarely, minima), but these do not include in- 
dividual measurements. These data can give a reasonable estimate 
for the period if sufficient successive maxima are available. 



2.2 Light curve analysis 

Mira light curves are often analyzed using the so-called O-C tech- 
nique (where O stands for the observed date of maximum and C 
for the calculated date). The O-C technique must be used with care 
when searching for secular period changes because it can be af- 
fected by period jitter: a sm all jitter can lead to larg e phase differ- 
ences over long time scales (Lombard & Koen 1993). In this paper 
we prefer the use of wavelet transforms. 

Wavelets are useful when the pulsation properties change over 
time, and have beeji^ise^l^ojjtudvjtjrig-^ 

mary et al. 1996; |Bedding et al. 1998J Kiss et al. 1999| ). We use the 
weighted wavelet Z-transform (WWZT jFoster 1996 ) developed at 
AAVSO specifically for uneven sampled data. 

We experimented with different values for the parameter c, 
which defi nes the trade off between time resolution and frequency 
resolution (Foster 1996), and settled on c = 0.005 as a good com- 



promise. More details of the appli cation of the WWZ tr ansform to 
long-period variables are given by Bedding et al. (1998). 



3 THE PERIOD EVOLUTION OF R HYA 
3.1 AD 1850-2001 

Fairly complete cov erage is available since 185 0, as dates of max- 
ima for 1850-1900 (|Cannon & Pickering 1909] M tiller & Hartwig 
1918), and as individual magnitude estimates since 1900. To al- 
low wavelet analysis, the pre-1900 dates of maxima and minima 
were arbitrarily assigned magnitudes of 5.0 and 10.0, respectively. 
In many cases only dates of maxima were available, and results 
from the wavelet analysis could only be obtained by inserting the 
missing minima. This was only done when it was clear that con- 
secutive maxima had been measured, in which case the date of the 
minimum was taken as being midway between the maxima. 

Fig. [| shows the wavelet plot for R Hya. The lightcurve is 
shown in the top panel, and the arbitrary magnitudes assigned to 
pre-1900 dates of maxima and minima are obvious. The second 
panel shows the WWZ transform, with the grey scale indicating 
the significance of each frequency as a function of time (see B ed- 
ding et al. 1998). Only a small range of frequencies is shown - 
there was no evidence for significant power outside this range. The 
third and fourth panels show, for each time bin, the semi-amplitude 
(in magnitudes) and period (in days) corresponding to the peak of 
the WWZ in the second panel. Note that semi-amplitudes are not 
available from these data prior to 1900. 

The period evolution in R Hya is clearly visible, with an over- 
all decline between 1850 and 1950, from 455 days to 385 days. 
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Figure 1. The light curve, using data taken from the AAVSO (1 1347 points), AFOEV (1489 points), BAAVS (4369 points) and VSOLJ (1509 points), binned 
to 10-day averages 



There is significant period jitter in addition to the decline. Before 
1900, some of the jitter may be due to the uncertainty in the dates 
of maxima. (Different determinations of the same maximum can 
typically differ by a few days to a week, but occasionally much 
more.) 

As made clear in Fig. ^, the period of R Hya is no longer 
decreasing ( 3reaves 2000| ). With hindsight, we can say that the pe- 
riod stabilized at ~385 days in about 1950, since which time the 
period jitter has been limited to the range 380-395 days. This jit- 
ter isjvidiirHhejTorma^ange_^ & Lombard 
1995). The period stabilization was preceded by a short phase of 
rapid decline, almost 10% within two decades. 

Another result is the behaviour in the semi-amplitude. There 
was a decrease from 2.2 mag to 1.7 mag between 1910 and 1950, 
closely mimicking the period evolution. The rapid period decline 
around 1940 is especially well matched by the amplitude, as is the 
constancy of the amplitude since 1950. Fig. | shows the close re- 
lation between the semi-amplitude and period. We have reported 
similar behaviour in the Mira variables R Aql, BH Cru and S Ori 
( Bedding et al. 200C ) and we speculated that, at least in some cases, 
the amplitude changes might cause the period changes via non- 
linear effects. For R Hya, the semi-amplitude reached a minimum 
around 1975, and has slowly increased again since. 



3.2 AD 1784-1850 

Table [l] lists all recorded dates of maxima before 1850, based 
on in formation given by tviuller & Hartwig (1918|), A rgelander 
(1869) and |Cannon & Pickering (1909[ ). The data are too patchy 
for wavelet analysis because the majority of maxima were not ob- 
served. Instead, we determined the period from observations of 
maxima that were consecutive or separated by only a few cycles. 
Since the period at 1850 is clearly established as ~450 days (see 
above), we start from that date and work backwards. 

Argelander's extensive observations gave maxima in 1843 and 
1848, consistent with a period close to 450 days and implying that 
three intervening maxima were missed. Olbers (who discovered the 
period evolution) observed maxima in 1818 and 1823, which give a 
period of about 460 days (assuming three missing maxima), which 
also matches the maximum observed by Schwerd in 1827 (assum- 
ing two missing maxima). Maxima observed by Piazzi in 1805 and 
1809 imply a period of 477 days, assuming two missing maxima. 
Finally, the maxima observed by Pigott in 1784 and 1785 were 485 
days apart and were presumably consecutive. 

Note that many of the unobserved maxima mentioned above 
would have occurred at times of the year when R Hya was not read- 
ily observable from Europe, as shown in Fig. ra. Together, these re- 
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Figure 2. The wavelet analysis for R Hya, 1850-2000. Shown are: the light curve, the frequency, the semi-amplitude of the main frequency component and 
its period. 

suits suggest that the period of R Hya was decreasing during 1784- 3.3 AD 1662-1712 
1850 at roughly the same rate as the post- 1850 decline, as shown 

in Fig. ^. Unfortunately, there is an 80-year gap in observations of R Hya 

prior to those by Pigott in 1784. The sparse observing record 
reflects the poor observability: at the star's southern declination 
(— 23deg), evening observations from Europe are feasible only in 
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Table 1. Pre- 1850 observations of R Hya. D ata taken from Miiller & 
Hartwig (1918) and Cannon & Pickering (190^ ) 
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Figure 3. Relation between period and semi-amplitude, from post- 1910 ob- 
servations. The trend visible in the upper right part of the plot, of decreasing 
amplitude with decreasing period, relates to the 1910-1950 data. The group 
of points in the lower part which show scatter contain the post- 1950 data. 



March-May (an early morning observation is required in winter), 
and even from Paris the star never reaches an airmass less than 3. 
But the large gap in the data also coincides with the depths of 
the Little Ice Age , with indication s for increased cloud cover over 
Northern Europe (Neuberger 197C). 

We are therefore left with the observations of R Hya made in 
the first decades after its discovery, which we now describe in (for- 
ward) chronologic a l order. The details are t aken from papers b y 
Argelander (1869| ), [vluller & Hartwig (1918| ) and ^offleit (1997| ). 
The first of these, in particular, contains a wealth of historical in- 
formation on several Mira variables. 

The first recorded observation of R Hya was by Johannes 
Hevelke (1611-16 89; latinized H evelius), who included it in his 
second catalogue (Hevelius 1679) but did not note any variability. 



The observations were made from Gdansk, Poland, on the evenings 
of Tuesday 18 and Wednesday 19 April 16620. 

Th e magnitude foun d by Hevelius is not certain. In his cat- 
alogue (Hevelius 1690b), he gives it as a 6th mag star, but ac- 
cording to_Mjiraldi_(see below) Hevelius observed it at 5th mag- 
nitude (Argelander 1869). Argelander states that he does not know 



how Maraldi obtained this value. However, Maraldi is likely to be 
correct^J^I^cWisjnd^aUai^^ 

1690a) (published shortly after this death) as of similar brightness 
to ip Hydra (my = 4.97). The included stars are consistent with 
a brightness limit i n this southerly region at my = 5. The chart 
of Hevelius (1690a I is reproduced in Fig. Q In the region below R 
Hya, a group of stars with magnitudes of 5.5 and fainter are lacking. 



J The lutheran Hevelius used the Gregorian calendar, made clear because 
he gives the days of the week of the observations. At this time, the Ju- 
lian calendar was still in use in Protestant parts of Europe, but Poland had 
adopted the Gregorian calendar in AD 1584, while neighbouring Prussia 
had done so in AD 1600. 
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a These are dates of observations rather than dates of maxima. However, 
given the magnitude range these observations could only have been made 
within 1-2 months of maximum. 

b There are two possible dates for this observations (see text). 
c Chandler (Astronomische Nachrichten 2463) derived the five maxima 
based on Maraldi's data. Dates in the table are as given by Argelander: 
Maraldi gives 1704 Marc h 14 and 1708 May 22. d The uncertain dates of 



maximum are as given by Miiller & Hartwig (1918). The uncertainties are 
due to the fact that the maximum was not covered, or in one case (1712) was 
difficult to observe because of the Full Moon. We estimate the uncertainties 
as 1—2 months. 

e Pigott mentions simultaneous observations by Goodricke which have not 
been published. 

f He also reports observations between 14 March and 4 May 1817, post- 
maximum, and between 13 Feb and 15 May 1822, also with an earlier max- 
imum. 

9 The date of this maximum is given by Schmidt (independent observation) 
as 3 (or 5) May. 



Cannon & Pickering ( 1 909t> gave the dis covery observation as the 
date of maximum, but Argelander ( 1 869 ) argued that the real max- 
imum occurred up to 2 months earlier or later. But if the star was a 
magnitude brighter than assumed by Argelander, it could have been 
closer to maximum^. 

Given this brightness limit, the variable U Hydrae should also 
be within the range of the catalogue. This carbon star varies be- 
tween visual magnitudes of 4.7 and 5.7. It indeed appears to be 
present in the c hart. There may be other, even older Chinese records 
of this variable ( frrtoffleit 1997 ), although its semi-regular variability 
was not discovered until 1871. We have not investigated U Hya fur- 
ther, but note that the stars shown in this region are also consistent 
with a limit of little fainter than 5.0. 

Geminiano Montanari, independently and before Hevelius' 
catalogue was published, observed R Hya while working at the 
Paris Observatory. While comparing Bayer's Uranometry with the 
sky, he noticed an unmarked 4 th magnitude star along the line con- 
necting if) and 7 Hya. Montanari did not publish the discovery and 
it is not known whether he noticed the variability. (Montanari had 
discovered the variability of j3 Per a few years earlier.) He entered 



4 Red stars can appear brighter to the naked eye in conditions of bright 
moon light, but the discovery date coincided with new Moon. 




the star with its magnitude on the map of Bayer. The date of his 
observation was 15 April in either 1670 or 1672 (see below). 

The variability of R Hya was first established by Giacomo Fil- 
ippo Maraldi (the nephew of Cassini). In 1702 he tried to re-identify 
R Hya based on Montanari's chart, but failed. But in March 1704 he 
observed the star and followed its appearances and disappearances 
until 1712 and identified maxima in 1704 and 1708rl According to 



Miiller & Hartwig (1918), Maraldi suspected a period of two years 
but, as they point out, this contradicts his own observations. They 
quote Pigott as deriving a period of 494 days in 1786 from his own 
and Maraldi's observations. As we can see in Fig. ^[ this period fits 
the five maxima of 1704-1712 very well. The period also agrees 
with Maraldi's failure to detect the star in 1702, when it would have 
been near minimum. 

The accuracy of the dates of maximum should not be over- 
stated. The high airmass worsens th e effects of the colours of the 
comparison stars. Argelander (1869) classifies the maxima accord- 
ing to accuracy, but for even the best determinations (1784 & 1785, 
1823, 1848) he estimates the uncertainty as 6-7 days. 

The evidence seems convincing that the period of R Hya dur- 
ing the time of Maraldi was about 495 days. As can be seen in 



Fig. q this indicates that the rate of period decrease was much less 
in the 18th century than in the 19th. 

We now turn to the two pre- 1700 observations, by Hevelius 
(in 1662) and Montanari (in either 1670 or 1672). The uncertainty 
in the year of Montanari's observation is unfortunate. The obser- 
vation was published by Maraldi, along with his own observations, 
in the Memoires de Paris (pour l'an 1706 and 1709), where the 
date was given twice as April 1672 and twice as 1670. However, in 
his own calculations Maraldi consistently used 1670. On the other 
hand, Montanari himself did not mention the star in a short paper of 
an academic speech from 1671 or 1672, describing several novae, 
which could fa vour the later datef Krgelander (1869|) and Mtiller & 
Hartwig (1918) preferred 1670, while |Cannon & Pickering (1909 ) 
gave 1672 as the more likely. 

Assuming that the observation by Hevelius was made close 
to maximum, any period close to that from Maraldi's data predicts 
a minimum in or around April 1672. On the other hand, accept- 
ing 1670 as the correct date and using a period of 496 days, we 
find an excellent fit with the observation by Hevelius and also with 
those of Maraldi (see Fig. ||). On their own, the two observations 
on 1662 and 1670/2 are consistent with an almost unconstrained 



The later observations by Maraldi were reported by Cassini (1740) 



6 The paper refers to a book in preparation on the 'Instabilita de firma- 
mento' but this never appeared. 
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Figure 5. The dates of observations are shown with a sine curve to indicate 
the proposed variability of R Hya. The plotted sine curve has a constant 
period of 496 days until 1770, declines at 0.4 days per year until 1800 after 
which the decline is 0.78 days per year. (This is not a unique fit, due to 
the paucity of data.) The solid curves indicate the parts of the light curve 
where observations before midnight would have been possible from Europe: 
March-June with m < 7. The dashed line extends this to include early 
morning observations. 



range of periods, and they can be made to fit any changing period. 
But the fact that Maraldi's confirmed period also fits these older 
observations leads to the plausible hypothesis that the period at the 
time of discovery was const ant, at about 496 days . We do not find 
support for the statement by Wood & Zarro (1981) that "four very 
old (1662-1708) and valuable dates of maximum ... show that the 
period was increasing." 



With this assumption of constant period, all observed maxima 
from 1662 to 1712 can be fitted to within 3 weeks (with the ex- 
ception of the poorly determined maximum in Nov. 1709, which is 
predicted 2 months earlier). Maxima would have occurred in Febru- 
ary 1662 and in March 1670, in good agreement with actual mea- 
surements. The lack of repeat observations by Montanari is also 
explained: the figure shows that the star would have been difficult 
to find for 2-3 years after his observation. 

3.4 The period evolution 

In summary, the period was approximately 495 days around AD 
1700, declining to 480 days by AD 1800, 450 days by AD 1850, 
420 days in AD 1900 and 380 days in AD 1950. The decline was 
almost linear, at 0.58 days/year: extrapolation suggests that the de- 
cline may have started around 1770, but it is also possible that the 
decline was initially slower and began earlier. Unfortunately, the 



decline probably began during the long gap in observations. The 
phases of constant period after and (possibly) before the decline 
suggest the possibility that the star has evolved from one (quasi) 
stable period to another. 

Fig. ^| shows how the observed dates of maxima fit with the 
period evolution. It is difficult to fit all observations with a purely 
linear period decline with a sudden onset. The fit used in the figure 
assumes a constant period until 1770, declining by 0.4 days per year 
between 1770 and 1810, with the decline increasing to 0.6 days per 
year after 1810. All dates of maxima before 1850 can be fitted well 
with this evolution. However, the constraints are relatively poor and 
equally good fits are feasible without assuming a gradual start to the 
period decline. Instead the result can be taken as evidence for some 
period jitter. The post-maximum observations of Olbers in 1817 
and 1822 are indicated as maxima at 01 Feb of those years. The 
figure shows that in both cases, the window of observability indeed 
fell post-maximum. 

The full period evolution is shown in Fig. ^. The dashed line 
is the fit proposed by Chandler (18961 in the third catalogue of 
variable stars. The sinusoidal component is not confirmed but the 
slope of his fit gives a good approximation until the period decline 
ended around 1950. 

A linear decline implies a constant rate of change in fraction 
per cycle, 5 = 1.6 X 10~ 3 . The time scale of the decline, defined 
as r = P/5, where P is the initial (longest) period, is r = 880 yr. 
This is an average time scale: the period evolution also showed 
significant jitter, with a fastest time scale (around 1940) of r m 
200 yr. 

Amplitude data, available since about 1890, show that the de- 
cline in period was accompanied by a decline in semi-amplitude, 
from 2.2 to 1.7 mag since 1905, or 5 mmag/cycle. The time scale 
for the amplitude decline, extrapolating back to 1770, is about 800 
yr, the same as for the period decline. The relation between period 
and amplitude is roughly linear (Fig. ^). The visual amplitude of 
an oxygen-rich Mira depends on the temperature variation during 
the pulsation, leading to the formation of molecules (TiO, VO) dur- 
ing minimum which strongly absorb at optical wavelengths (Reid 
& Goldston 2002). A relation between amplitude and period could 
therefore be strongly non-linear, but this is not seen in R Hya. 



4 STELLAR PARAMETERS 

The main uncertainty in the mass and luminosity of R Hya de- 
rives from its uncertain distance. Unfortunately, the Hipparcos par- 
allax is a non-detection: 1.6 ± 2.4 mas. W hitelock, Marang & Feast 
(2000) found a distance of 140 pc, derived by placing R Hya on 
the Mira period-luminosity (P-L) relation. Eggen ( 1 985 , 1966) 



reported a proper-m otion companion which g ave a distance modu- 
lus of 6.1 (165 pc). lura & Kleinmann (1992) favoured a distance 



of 1 lOpc, based on a P-L relation. Only the value of Eggen is con- 
sistent with Hipparcos at the 2-a level: these two are also the only 
direct measurements. 

Using d — 165 pc, the luminosity is L — 1.16 x 10 4 L®. The 
mean of the Eggen and Hipparcos distances, 400 pc, would already 
yield a luminosity above the classical AGB limit (Mbol = —7.2: 
this limit can be exceeded in th e case of hot bottom burning, bu t 
only for very large core masses: Bloecker & Schoenberner 1991 ). 
For this reason we will use Eggen's distance in the discussion to 
follow. 

Eggen (1985) argued that R Hya is located within the Hyades 
supergroup, with an age of 5-10 x 10 s yr. This would imply a pro- 
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Figure 6. The period evolution of R Hya between 1662 and 2001. The first point is uncertain; the p eriod is well dete rmined from 1704 onwards. Extrapolation 
of the linear decline suggests the decline began around 1770. The dashed line is the fit proposed by ^handler (189^ ). 



genitor mass for R Hya around 2 M© . The presence of technetium 



( ja Tc) in R Hya (Little et al. 1987) shows that the star is in the 



thermal pulsing phase of the AGB (e.g., Lebzelter & Hron 1999) 



this element is dredged-up during the thermal pulses but has a half 
life of ~ 10 5 yr, several times the interpulse time scale. Its abun- 
dance slowly increases during the TP-AGB. Tc is found in 15% of 
semiregulars but 75% of long-period Miras. 

Iru^are^^hotonrtr^^as_^ejjorte^^v_Whitelock et al. 

(2000): phase-averaged magnitudes are (J, H, K, L)o = 
-1.07,-2.05,-2.45,-2.88) and the bolometric magnitude 
is mboi = 0.64. The infrare d colours indicate an effective 
temperature of T cff = 2830 K ( feast 1996]). The V - K ^9.5 
is consisten£jvid^his^e^rr£erat^ 

et al. (1995) obtain a lower T cB = 2680 ± 70 K by fitting to the 
flux distribution between 1.04 and 3.45 microns. They also derive 
T c ft from an angular diameter measurement: they found 2570 K or 
2760 K, assuming fundamental pulsation mode or first overtone, 
respectively. 



The luminosity, with the temperature derived by Feast (1996), 
yields a radius of R ~ 450 R© . For our adopted distance, this pre- 
dicts an angular diameter of 26 mas. Th e angular diameter has been 
measured at 902 nm as 34. 1 ± 3.4 mas ( Haniff et al. 1995 ), assum- 
ing an uniform disk, yielding a large radius of R — 590 R© . A cor- 
rection for limb darkening and molecular opacities brings the value 
down to about 500 R© for first overtone models; for the fundamen- 
tal mode the effect is much less. Recently, Tuthill (priv. comm.) 
measured a near-infrared diameter of 24 mas, in much better agree- 
ment with the prediction above. Part of the difference between the 



two observations may be due to photospheric extensions which can 
be si gnificant at 9 02 nm: the if -band is likely to be less affected by 
this (Feast 1996). In addition, the earlier observation took place 



close to minimum (Tuthill, priv. comm), which in Miras occurs 
when the star is largest. 



Whttelock et al. (2000) and Feast (1990 assumed that R Hya 



is presently located on the Mira P—L relation. However, this re- 
quires a distance (140 pc), outside the 2-cr confidence limit of Hip- 
parcos. Given its period history, a location on the narrow P~L rela- 
tion may not be expected. The distance assumed in this paper would 
put R Hya slightly above or to the left of the relation, perhap s be- 
tween the Mira and SR branches (Bedding & Zijlstra 1998). The 
Hipparcos distance places the star significantly above the relation, 
a location in common with O-rich LM C Miras with P > 420 days 
(jFeast et al. 1989] fcijlstra et al. 1996]). 



5 THE PULSATION 

The gradual change in the period of R Hya implies that its pulsation 
mode has remained constant; its evolution is therefore related to a 
change in the stellar parameters. 

The pulsation mode of R Hya is an open question, as it is for 
all Mira variab les (|Wood 1990|; IvVhitelock & Feast 2000|; Y a' Ari & 
Tuchman 1999). Neither is it proven that R Hya exhibits the same 
pulsation mode as other Miras. The radius derived above is consis- 
tent with either the fundamental mode or with first overtone (e.g., 
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Whitelock & Feast 2000 : R Hya falls in between the two modes in 



their Fig. 1). 

The pulsation equation, which relates the period P (in days) 
to the radius R and mass M (in solar units) is given by: 



for first overtone pu lsators, where the pulsation constant Q ~ 0.04 
(|Fox&Wood 1982b, or 



for fundamental mode pulsators (Wood 199C). These equations 
yield masses for R Hya of 0.74 Mq and 3.0 M , respectively. The 
large mass required for the fundamental mode provides an argu- 
ment for the first overtone, or alternatively for questioning whether 
the measured angular diameter is identical to the pulsational diam- 
eter. 

The two equations both imply that the period evolution was ac- 
companied by a change in radius: the radius would have decreased 
by 14-18%, depending on pulsation mode. The pre-1770 radius 
would have been about 520 Rq . 

There are no direct observations to show how T c ff and L 



changed during the period evolution. Wood & Zarro (1981) fitted 



a luminosity decline of 20% , based on the assumption that R Hya 
underwent a thermal pulse. Ya'Ari & Tuchman ( 199€ ) presented 



a different model for period evolution (see below) which does not 
require a change in luminosity. The lack of information on the lu- 
minosity evolution does not allow us to test these two models. 
The P-L relation derived from LMC Miras is given by: 



Mboi = -3.001ogP + a 



(Feast 1996). This predicts a luminosity decrease of 25% for the 



log P = 1.5 log R - 0.5 log M + log Q, 



(1) 



log P = 1.949 log R - 0.9 log M - 2.07 



(2) 



(3) 



period decline of R Hya. However, this should be taken as an up- 
per limit, as R Hya is unlikely to have evolved along this relation: 
the P-L relation is not an evolutionary sequence but rather a se- 
quence of stars with diffe rent progenitor masses and metallicities. 
The evolutionary tracks of Vassiliadis & Wood (1993| ) cross the P— 
L relation at almost constant luminosity, while the Whitelo ck evo- 
lutionary track found in globular clusters ( Whitelock 1986 ) is also 
shallower. But short-term evolution, such as that shown by R Hya, 
may not follow these sequences either. 

Combinin g the relation between colour and period of W hite- 
lock et al. (2000), 



(J - K) = -(0.39 ± 0.15) + (0.71 ± 0.06) logP, 



with the Tog-colour relation from Feast (199€ ) 



logTeft = -0.59(J - K)o + 4.194, 



(4) 



(5) 



yields an increase of T e g for R Hya of 10%, i.e., from 2570 to the 
present 2830 K. Combining this with the radius change from the 
pulsation equation gives the counterintuitive result that the lumi- 
nosity of R Hya has increased by 5% rather than decreased. Given 
the slope of the P-L relation, this suggests that the slope of the 
temperature calibration used here is too steep. The temperature cal- 
ibration averages Mira and non-Mira M-type stars. Using only Mira 
variables gives a more shallow relation: 



logTeff = -0.474(J - K) + 4.059, 



(6) 



which gives a 3% decrease in luminosity. These relations suggest 
a negligible change in luminosity. The assumption that R Hya re- 
mained on the AGB colour relations (at constant L) may be more 
realistic than that of R Hya following P-L relations, which predict 
decreasing L. 
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6 DISCUSSION 
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6.1 Mass loss evolution: winds of change 

There are strong observational relations betwe en stellar parameters 
and mass loss on the TP-AGB. |Blocker (1995 ) proposed a variation 
of the Reimers mass-loss equation: 



M B = 4.83 x 10~ 13 M' 21 L 2 7 



( — ) 

KM J 



(10) 



whe re the last term comes from the Reimers equation ( Reimers 



1975). Vassiliadis & Wood (19931 used a very different formula- 
tion: 



log Afvw = -11.4 + 0.0123P 



(11) 



for stellar winds below the radiation momentum limit. Both rela- 
tions predict a change in M for R Hya during its recent evolution. 
The Blocker equation predicts, for a change in radius of 15% and 
in luminosity of 20%, that the mass-loss rate would have declined 
by a factor of 3. The decline is governed mainly by the luminosity, 
for which we have used the most extreme estimate. If the luminos- 
ity has remained consta nt, the mass-loss decline w ould be much 
smaller. In contrast, 



the 



Vassiliadis & Wood (1993) relation pre 



diets a much steeper decline, by a factor of 20 independent of any 
luminosity evolution. Their relation also predicts a decline of the 
wind expansion velocity from 14 to 8 km/s. (Both relations are used 
to model evolutionary tracks and may not describe the short-term 

changes in R Hya.) 

Hashimoto et al. (1998) drew attention to the peculiar IRAS 
spectrum of R Hya, which shows a dust continuum without sil- 
icate feature (class In). Silicate emis sion forms close to the star 
and its lack indicates a detached shell. Hashimoto et al. (1998) de- 



rived an inner radius of 60i?», based on a distance of HOpc and 
R* — 700 R.0 . To first order, the inner radius scales with luminos- 
ity. Scaling to Eggen's distance gives Rj = 6 x 10 15 cm. For an 



expansion velocity of 7.5 km/s ( Wannier & Sahai 198Q , Ri corre- 
sponds to 250 yr BI (before IRAS). This would put the decrease of 
the mass-loss rate around AD 1750. 

The uncertainty in this calculation is significant, not least be- 
cause the fit assumes a sudden end to the mass loss, while a gradual 
decrease is more likely. (The mass loss has not ceased completel y, 
as shown by the presence of an SiO maser (Snyder & Buhl 1975).) 



The outer radius indicated by the fit is 10 cm, although uncertain. 
This corresponds to an age of 3000 yr. 

The (pre-1770) mass-loss rate derived by Hashimoto et al. 



M 



, which is 



(1998), scaled to d = 160 pc is M w 3 X 10" 
low for a long-period Mira (compare Afvw = 5 x 10 -6 Mq yr -1 ). 
A value around 10~ 7 Ma yr" 1 is obtained from the CO(2-l) mea- 



surements (Wannier & Sahai 1980- Hashimoto et al. (1998) did 



not give limits on the present-day mass loss, but the lack of any 
silicate suggests that the decline was more than predicted for Mb, 
and perhaps closer to the prediction for Mvw- 

To esti mate the required decrea se in M, we have repeated the 
model fit of Hashimoto et al. (1998). With a single wind, we con- 



firm the mass-loss rate and cavity size required by the LRS spec- 
trum. If we fill the cavity with a lower-density wind, a weak sili- 
cate feature re-appears. Only with the new wind at least 10 times 
less dense can we fit the spectrum. This is a much larger decrease 
than predicte d by Blocker's formalism b ut is in agreement with the 
prediction of Vassiliadis & Wood (19 93). The strong decrease pre- 
dicted by Vassiliadis & Wood (1992) appears to be confirmed for 
RHya. 

Interestingly, the IRAS 60-^im image shows a detached shell 
around a bright point source, with an inn er radius of 1-2 arcmin 
(1.5-3 10 17 cm). Hashimoto et al. (1998 ) argued that this gap is 



inconsistent with their model, with an inner radius that is far too 
large, and they cautioned that the deconvolution procedure used 
(Pyramid Maximum Entropy) can give artifacts in the presence of 
a bright point source. However, the possibility should be considered 
that this ring represents a much older mass-loss event. Its inner ra- 
dius indicates that this mass-loss phase was interrupted ~ 5000 yr 
ago. 



6.2 Real-time evolution 

For the time scale on which R Hya evolves, two models have been 
described in the literature that fit its period evolution. 



6.2.1 Post-thermal-pulse evolution 

A thermal pulse occurs when sufficient helium has built up from the 
ashes underneath the hydrogen burning layer. The TP gives a strong 
modulation of the stellar luminosity. At first, the luminosity spikes 
over a time scale of ~ 10-100, yr. Then the luminosity reaches a 
short-liv ed plateau at a level above the h ydrogen burning luminos- 
ity (e.g., Boothroyd & Sackmann 1985), followed by a decline on 



a time scale of a few hundred years. The luminosity continues to 
drop slowly during quiescent helium burning, reaching around 1/3 
of the hydrogen burning luminosity. This phase lasts about 10% of 
the TP cycle. Finally, helium burning ceases and the hydrogen layer 
re-ignites, quickly recovering the pre-pulse luminosity. The period, 
and also the mass-los s rate, mimic the luminosity evolution ( Vas- 
siliadis & Wood 1993; plocker 1995j ). Roughly speaking, the TP 
phase lasts 10 2 -10 3 yr, the helium-burning phase 10 3 -10 4 yr and 
the quiescent hydrogen-burning phase 10 4 -10 5 yr. Detached shells 
aro und AGB stars are c ommonly interpreted in terms of the TP cy- 
cle (jZijlstra et al. 1992;). 

Wood & Zarro (1981) located R Hya within the earliest post- 
TP evolution, when the luminosity shows the steepest drop. In their 
fit, the peak luminosity would have occurred around 1750 and the 
period (and luminosity) during the Hevelius-Montanari-Maraldi 
observations would have been increasing. We have shown that there 
is no evidence for a period increase, although it cannot be ruled 
out either. Sadly, there are no observations during the crucial phase 
around 1750. A near-constant period during 1662-1784 could still 
be accommodated in their mo del by assuming the p ulse occurred 
50 yr earlier than assumed by Wood & Zarro (1981), placing the 
peak luminosity plateau around 1700. Their model also predicts a 
slowing of the luminosity decline around the present time, which is 
consistent with the observed lack of evolution since 1950. 

The TP model fits the time scale and period decline well. A 
concern is that it places R Hya within a unique 100-500 yr phase of 
the TP cycle, corresponding to only 1% of the cycle. The likelihood 
of this occurringjn the 3 rd brightest Mira on the sky is small. Sterne 



skyj 



et al. (1937) found continuous period changes in 2 out of 377 well- 
studied Miras, which is in agreement with this TP-phase. (A few 
more Miras are now known with large period evolution: Bedding et 
al., in preparation). 

The duration of the high mass-loss phase pre-1770 may be 
more difficult to reconcile with the TP model. If this phase traced 
the peak of the pulse, a duration of ~ 10 2 yr would be expected, 
while if it traced the phase of quiescent H-burning it should have 
lasted ~ 10 4 yr or longe r. Both the model and the IRAS images of 
Hashimoto et al. (1998) suggest it lasted for several 10 3 yr, which 
is consistent with neither. 

The evidence for an earlier mass-loss interruption also raises 
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a problem. With a time difference of ~ 5000 yr, it is not possible 
to relate bo th to a thermal puls e. If the first event was due to a ther- 
mal pulse (Zijlstra et al. 1992), R Hya would presently be nearing 
the end of the helium-burning phase or have recently re-entered the 
higher-luminosity hydrogen-burning phase, a phase with a much 
slower luminospity evolution. For the TP-model, it would be im- 
portant to investigate whether the detached ring in the IRAS image 
is real or could be explained as an imaging artifact. 



6.2.2 Envelope relaxation 

Mira pulsations are intrinsically non-linear. The period may depend 
on the amplitude of the pulsation, affecting either the radius R or 
the pulsation constant Q. The fact that the amplitude and period of 
R Hya show e vidence for simultaneous evolution (see also M attei 
& Foster 20 00) could show the presence of such a non-linearity. 
Wood (1976) suggests that small variations in Mira period are most 



easily explained by an alteration in the envelope structure near r — 
0.8R,. 

The effect of non-linearity is studied by Ya'Ari & Tuchman 
(1996), who calculated the pulsational stability over a much larger 
number of cycles than had been done before. In their models, fol- 
lowing an induced perturbation, the star pulsates in the first over- 
tone for ~ 200 yr. During this time the growth rate of the funda- 
mental mode is small but non-zero. Once the fundamental mode 
begins to dominate, a re-arrangement of the envelope structure oc- 
curs, with entropy transported downward. The period of the funda- 
mental mode slightly declines when this mode first dominates, but 
during the change of the stellar structure the fundamental period 
declines over a period of ~ 150 yr. Their model closest to R Hya 
is model D, where the period first declines to 495 days, and during 
the restructuring declines to 330 days. This change is a little larger 
than seen in R Hya but occurs on a very similar time scale (but note 
that the model star has a much lower luminosity than R Hya). 

The strong points of the model are that the onset, time scale, 
and eventual stabilization can all be explained. However, it requires 
that the star is initially in a non-equilibrium state and the cause of 
this is open. The average luminosity is constant during the period 
evolution. 



6.2.3 Cause and effect 

In the TP model, there is a clear cause for the change in period: 
the declining luminosity causes a reduction in the stellar radius, 
which causes the period to become shorter. In the envelope relax- 
ation model, what triggers the mode switch is an open question. 
One possibility is the effect of weak chaos. Icke et al. (1992) 



showed that the outer layers of the star can lose track of the un- 
derlying pulsation and become trapped in 'islands of stability'. The 
effect is strongest for stars that have reduced envelope masses, and 
has been i nvoked to explain the mode switching in R Dor (B edding 
et al. 1998). In the model of |Icke et alj , there is an underlying piston 
moving with constant frequency. In real Miras, the non-linearity 
discussed by ya'Ari & Tuchman (1996| ) implies that if a star is 
caught in an island of enlarged radius, over time the inner structure 
of the star could be affected by this. This could act as the trigger 
for the mode evolution. 

Interaction betwee n the star and its extende d atmosphere may 
also have some effect: Hoefner & Dorfi (1997) have shown that 
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Figure 7. Relation between the separation of the rings observed in post- 
AGB stars, and the location of the rings, both in arcsec. The squares corre- 
sponds to poitions midway between the density enhancements. 



We find a clear relat ion between amplitude and period for R 
Hya. Bedding et al. (2000) have suggested that the change in ampi- 



tude may act as the cause of the period change. 



6.3 Rings 

Several post-AGB stars and one AGB sta r show concentric rings 
seen in reflected light (Kwok et al. 2001). The separation of the 
rings (or arcs: only the illuminated parts are seen) correspond to 
time scales of about 500 yr. The thicknesses of the rings correspond 
to 0.1-0.5 of the separation, and the density enhancement in the 
rings is at least 30%, but could be larger. The obvious explanation 
of these rin gs is that the mass -loss rate showed a fluctuation on this 
time scale (Sahai et al. 1998). However, the only effect known to 



feedback from atmosphere on the star can affect the cycle-to-cycle 
amplitudes. 



modify the AGB star properties on this time scale is the thermal 
pulse, and this could only lead to a single ring. 

The time scale for the period decline in R Hya is remarkably 
similar to the time scale of the rings. The strong decline in mass- 
loss rate following the onset of the period decline makes it the only 
observed Mira behaviour which can explain the rings. However, 
this requires the evolution discussed in this paper to be periodic. 
The fact that the period has now stabilized allows for the possibility 
that it will at some time increase again, but there is at present no 
direct evidence that the period evolution is periodic. 

Of the models for the R Hya evolution discussed above, only 
the relaxation model combined with a periodic or stochastic trigger 
could lead to the formation of multiple rings. The observed separa- 
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tion in the rings is not fully regula r: this is shown in Fig. (7| using 
data taken from Kwok et al. (2001 ). The separation can vary by as 



much as a factor of 2 (although in a few cases an intervening ring 
may have been missed). There is also a clear indication for an in- 
crease of rings separation with distance from the star. This implies 
that the event causing the rings occured at decreasing time intervals 
as the star approached the end of the AGB, 

The chaotic behaviour predicted by [eke et al. (199^ ) increases 
as the envelope mass reduces, and this behaviour fits both the irreg- 
ularity and the increasing frequency of the mass-loss epsiodes. But 
its effect on the mass loss is not clear, and it is not proven (although 
possible) that this chaotic behaviour can act as a trigger for an R 
Hya-type event. 

The TP model makes a very clear prediction for the future pe- 
riod evolution. Further monitoring of R Hya is therefore important: 
a continuing but slow decline would agree with the TP model. If, 
in contrast, the period is found to increase again, this would rule 
out the TP model and make a connection with the post- AGB rings 
more likely. 



7 CONCLUSIONS 

We have studied the period evolution of R Hya, using both mag- 
nitude estimates for the light curve and old data giving dates of 
minimum/maximum. The wavelets are shown to be a powerful tool 
to analyse such datasets. The main results are 

(i) The period of R Hya has declined continuously from 495 
days to 385 days, between approximately AD 1770 and AD 1950. 
Before 1770 there is no evidence for period evolution, while after 
1950 the period has been stable, showing at most minor period jit- 
ter. The evolution gives the impression of a change between two 
relatively stable configurations. We do not confirm the suggestion 
that prior to 1770 the period was increasing. 

(ii) The amplitute (available after 1900) closely followed the pe- 
riod evolution, declining at first but becoming stable after 1950. A 
relation between amplitude and period is typical for a non-linear 
pulsation. 

(iii) The most likely distance is 165 pc, giving a luminosity of 
1.16 x 10 4 L0. The likely progenitor mass is around 2Mq. The 
star is located on the thermal-pulsing tip of the AGB. 

(iv) The period change indicates a decrease in stellar radius. 
The luminosity and temperature change is less secure. Assum- 
ing the star remained on the fiducial AGB relations, the tempera- 
ture change may have been 10-20%. Various luminosity-dependent 
AGB relations predict changes in the luminosity ranging from 25% 
decrease to 3% increase. Given the uncertainty whether R Hya sat- 
isfied such relations during its period decline, and the fact that dif- 
ferent relations do not even agree on the sign, it is not possible to 
confirm that the luminosity decreased: a constant luminosity is a 
significant possibility. 

(v) The IRAS spectrum shows that mass-loss rate has recently 
declined by a factor of at least 10. A model of the IRAS spectrum 
shows that the mass-loss decline occured about 250 yr ago. This is 
in good agreement with the onset of the period decrease and sug- 
gests the two effects are correlated. The pre- 1770 mass-loss rate 
was M«3x 10~ 7 M yr -1 . A large detached IRAS shell sug- 
gests an earlier phase of high mass loss, ending about 5000 yr ago. 
The post- 1770 decline agrees with the M—P relation of Vassiliadis 
& Wood (1992) but is much larger than predicted from the mass- 
loss formalism of Blocker (1995). 



(vi) Two models can explain the behaviour of R Hya. First, a 
recent thermal pulse, occuring shortly before the discovery. This 
also can fit the constant period since 1950. The second model is 
envelope relaxation, where the non-linearity of the Mira pulsation 
causes a change in the entropy structure of the star. Both the pe- 
riod evolution between two semi-stable states and the time scale of 
the change are reproduced. There is at present insufficient data to 
decide between the two models. 

(vii) The evidence for a strong effect on the mass loss raises 
the possiblity of a connection with the circumstellar rings observed 
around some post- AGB stars. The evolution seen in R Hya is the 
only observed effect in Miras which has the correct time scale. Icke 
et al. (1992) show that Mira period instability increases as the enve- 
lope mass decreases. This would place such events at the tip of the 
AGB, and would agree with the observations that the time scales 
between 'ring' events decreases with time. However, a mechanism 
to translate this chaotic envelope behaviour into structural (period) 
evolution of the Mira is lacking. 

(viii) Further monitoring of R Hya is recommended. The 
thermal-pulse model makes a strong prediction for its future pe- 
riod evolution. If, on the other hand, the period at some time would 
increase again, this would rule out this model and also make a con- 
nection with the post- AGB rings more likely. 

Changes in Mira properties were already known on a cycle-to- 
cycle basis, and on time scales of 10 4 yr, which is the thermal-pulse 
time scale. R Hya shows that significant evolution can also occur 
on intermediate time scales of order 10 2 -10 3 yr. 
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